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Abstract 

We discuss three-body recombination of ultra-cold atoms to a weakly bound 
s level. In this case, characterized by large and positive scattering length a 
for pair interaction, we find a repulsive effective potential for three-body colli- 
sions, which strongly reduces the recombination probability and makes simple 
Jastrow-like approaches absolutely inadequate. In the zero temperature limit 
we obtain a universal relation, independent of the detailed shape of the inter- 
action potential, for the (event) rate constant of three-body recombination: 
ttrec = 3.9fia^/m, where m is the atom mass. 
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Three-body recombination, the process in which two atoms form a bound state and a 
third one carries away the binding energy, is an important issue in the physics of uhra-cold 
gases. This process represents the initial stage in the formation of clusters intermediate in 
size between individual atoms and bulk matter. Three-body recombination limits achievable 
densities in high-field-seeking spin-polarized atomic hydrogen and in trapped alkali atom 
gases (see and references therein) and, hence, places limitations on the possibilities to 
observe Bose-Einstein condensation in these systems. 

Extensive theoretical studies of three-body recombination in ultra-cold hydrogen [jl|,|^ 
and alkalis showed that the rate constant of this process, Orec, strongly depends on the 
shape of the potential of interaction between atoms and on the energies of bound states in 
this potential. In alkalis the recombination is caused by elastic interatomic interaction, and 
in the zero temperature limit a^^c varies approximately as 0, where a is the scattering 
length for pair interaction. 

All these studies, except one in spin-polarized hydrogen (see 0), rely on Jastrow-like 
approximations for the initial-state wavefunction of three colliding atoms. Recent progress 
in the quantum three-body problem for the case where only zero orbital angular momenta 
of particle motion are important opens a possibility for rigorous calculations of three- 
body recombination in ultra-cold atomic gases. In this Letter we consider the extraordinary 
case of recombination (induced by elastic interaction between atoms) to a weakly bound s 
level. The term "weakly bound" means that the size / of the diatomic molecule in this state 
is much larger than the characteristic radius of interaction R^. (the phase shift for s-wave 
scattering comes from distances r ^ R,,). In this case the scattering length is positive and 
related to the binding energy Eq by (see, e.g., 0) 

a = lij ^msQ r-^l^ Rf. (1) 

(m is the atom mass), and elastic (s-wave) scattering in pair collisions is resonantly enhanced 
at collision energies E -^e^. As we show, large a and / imply a rather large recombination 
rate constant a^-^c- At the same time, for large positive a we find a repulsive effective 
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potential for three-body collisions, which strongly reduces a^ec- In the limit of ultra- low 
initial energies E <^ sq we obtain a universal relation independent of the detailed shape of 
the interaction potential: ctrec = 3.9ha^/m. 

The dependence a^ec oc can be understood from qualitative arguments. For atoms of 
equal mass the energy conservation law for the recombination process reads 

Sh'^kj/Am = eo, (2) 

where kf ~ 1/a is the final-state momentum of the third atom relative to the center of mass 
of the molecule. Recombination to a weakly bound s level occurs in a collision between two 
atoms, when a third atom is located inside a sphere of radius / ~ a around the colliding 
pair. For such locations of the third atom, characterized by a statistical weight w ~ nl^ {n 
is the gas density), this atom and one of the colliding atoms form the weakly bound state 
with probability of order unity. The number of recombination events per unit time and unit 
volume, i^rec = ctrec^^, cau bc estimated as ii?av{nl^), where a = Sna'^ is the cross section 
for pair collisions. One may put velocity v ~ hkf/m, which gives a^^c ~ Sirha'^/m. 

One can also understand qualitatively the existence of a repulsive effective potential for 
three-body collisions and the reduction of arec- In the mean field picture the interaction in 
a three-body system at (maximum of the three) interparticle separations r ^ Re can be 
written as Anfi^n^a/m^ where ra^, ~ is the "particle density" inside a sphere of radius 
r. For a > this interaction is repulsive, which makes the statistical weight w smaller than 
nl^ and decreases the numerical coefficient in the above estimate for ctrec- The tail of the 
three-body effective potential at r ^ a was found in P). Arguments clearly showing the 
absence of any "kinematic" repulsion independent of the value and sign of a are given in . 

A particular system that should exhibit three-body recombination to a weakly bound s 
level is a gas (or a beam) of helium atoms. The He-He potential of interaction V{r) has a 
well with a depth of 11 K. There is only one bound state in this well, with orbital angular 
momentum j = and binding energy Eq ^ 1.3 mK (see [^] and references therein). The 
scattering length a ~ 100 A found for this potential satisfies criterion (|l|). The existence 
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of the He2 dimer, the world's largest diatomic molecule (/ ~ 50 A), has been established 
experimentally Another system which is likely to have three-body recombination to a 
weakly bound s level is spin-polarized metastable triplet helium, a gas of helium atoms in 
the 2^5* state with spins aligned. The interaction potential |]10| for a pair of spin-polarized 



B.e{2^S) atoms supports an s level with binding energy ~ 2 mK, which leads to a ~ 100 A 



and important consequences for the decay kinetics of this system ||TT |. 

We confine ourselves to three-body recombination of identical atoms at collision energies 
-E ^ to a weakly bound molecular s level. In this case the recombination rate constant 
Qfrec can bc found from the equation 



27r rd^kf ,o./3^2H 



n 
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Here n^/Q stands for the number of triples in the gas, Tif = J iIjiViIj^P* d^xd?x' is the T- 
matrix element for three-body recombination, the coordinates (x, x') are specified in Fig.l, 
ipi is the true wavefunction of the initial state of the triple, and ip^p is the wavefunction of 
free motion of the third atom relative to the center of mass of the molecule formed in the 
recombination event. Both ipi and ifj^p can be written as a sum of three components, each 
expressed in terms of one of the three different sets of coordinates (see Fig.l): 

A = V^(x, y) + ^(x', y') + 7/;(x", y"), (4) 

= (l/v^)[0(x,y) + 0(x',y') + 0(x",y")], (5) 
0(x,y) = V^o(a;)exp(zk/y), 

where ipo is the wavefunction of the weakly bound molecular state. The interaction between 
colliding atoms is regarded as a sum of pair interactions V{r). The quantity V is the part of 
the interaction which is not involved in constructing the wavefunction i.e., if the molecule 
is formed by atoms 1 and 2 (the first term in Eq.^, then V = V{ri — r^) + V{r2 — r^), et 
cet. Using Eq.(||), 

Ti/= 2v^^^xrfVV^o(a;)cos|^^y (x')exp(-ikjx')V^i. (6) 



The initial wavefunction of the triple is best represented in hyperspherical coordinates. 
The hyperradius, defined as p = + 2?/^/3)^/^, is invariant with respect to the trans- 

formations x,y ^ x',y' — > x",y". The hyperangles are defined as a = arctan(v^a;/2?/), 
and similarly for a' and a". For E <^ sq only zero orbital angular momenta of the particle 
motion are important, and the wavefunction ip can be written as |^ 

^ ^ y-:^A(p)_$AKp)_ /^N 
^ ^/Q sin a cos a 

The functions $A(a,p) are determined by the equation 

|^+7^^(v^Psi^«)/''PA(«,p)+^^^^^_^l rfa'$,(a',p))=A(p)$,(a,p), (8) 

with boundary conditions $a(0,p) = $A(7r/2,p) = and normalization Jq'^^^ \^\{pL^ p)^da = 
7r/4. The sum in Eq.(0) is over all eigenvalues A corresponding to three free atoms at infinite 
interparticle separation. At ultra-low collision energies the lowest such A(p) alone gives a 
very good approximation, and we can confine ourselves to this A. Then the function -Fa(p) 
can be found from the (hyper)radial equation in which the quantity A(p) serves as an effective 
potential . Under the condition E <^ ai interparticle distances much smaller than their 
De Broglie wavelength this equation reads 

( Sa^AWziU.W.O. (9) 



y dp"^ pdp p^ 

The function F\{p) should be finite for p ^ and is normalized such that F\{p) 1 for 

p — s> OO. 

In our case the pair interaction potential V{r) supports a weakly bound s level, and the 
scattering length is positive and much larger than the characteristic radius of interaction 
for this potential. For p^ R^ the function (^x{a,p) takes the form (cf. 0) 

o(p)a [(v^p/a) sin(7rv^/2)xo(y2pa) + (8/ V3)sin(7rv^/6)1, a <i?e/p , , 
*A(a,p) = <{ ^ (10) 

g{p)sin -\/A (a — 7r/2) , a > -Re/p, 



where g{p) = [1 + sin(7rV^)/7r-\/A] ^^"^ and Xo('^) is the solution of the Schrodinger equation 
for the relative motion of a pair of particles. 



normalized such that Xo ~^ 1 — a/r asr oo. Matching the wavefunctions (|T0D at a = 
Re/p -C 1, to zero order in Re/p we obtain the following relation for A(p) at distances 
p»/?e (cf. i): 

For p ^ a this equation yields A(p) = 4 + 48a/-\/27rp, and thus the potential term in Eq.(^ 
varies as a/p^. Eq.(|T^ is universal in the sense that A depends only on the ratio p/a, but 
not on the detailed shape of V{r). The same statement holds for Fx{p) at distances p ^ Re- 

For infinite separation between particles, i.e., for p ^ oo and all hyperangles larger than 
Re/p, we have \/X ~ 2 and ^x{a,p) ~ sin 2a. Accordingly, from Eq.(^ with Fx{p) — > 1, 
each ijj in Eq.@ becomes equal to y^2/3, and the initial wavefunction ipi ^/6. 

The "effective potential" A(p) and the function F\{p) for three ground-state He atoms 
(a ^ lOOA) are presented in Fig. 2 and Fig. 3. The potential V{r) was taken from [§. 
For p ^ lOOA our numerically calculated A(p) coincides (within 10%) with that following 
from Eq.(|l^), ensuring a universal dependence of Fx on p/a. As A(p) is repulsive. Fx is 
strongly attenuated at p < a (see Fig. 3). This leads to a strong reduction of ipi when all 
three particles are within a sphere of radius ~ a. 

We first consider the theoretical limit of weak binding, where the scattering length a and 
the binding energy Eq are related by Eq. (|1|) , the wavefunction of the bound molecular state 
at distances x ^ R^ is 

11 / x\ 
= -i==-exp — , (13) 
V27rax \ a/ 

and the final momentum kf = 2/\/3a. From Eq.(|l3D one can see that the distance between 
the two atoms which will form the bound state should be of order a. To take away the 
binding energy the third atom should approach one of them to a distance of order Re. The 
main contribution to the integral in Eq.(^ comes from distances a; ~ a and a;' ~ i?e <^ o,- 



Therefore we may put p ~ y2/3x, a = a" ~ vr/S, and a' ~ \/?>x' /2x. Then the initial 
wavefunction takes the form 

V'i ~ (l/V3)xo(x')^A(V2a;/y3a), (14) 



with Fx{z) = zFx{z)g{z)?,m.{\J \{z)t: /2) and z = p/a. Putting kjx'^^O and using 
^S'x'V{x')xQ{x')=AT[h^a/m, from Eq.® we obtain Tif = AM^^'^Ti^a^/'^G/m, where 

G = dz 

The main contribution to this integral comes from z ~ 1 (p ~ a), where A and Fx (and, 
hence, Fx) are universal functions of p/a. Therefore G is a universal number independent of 
the potential V{r). Direct calculation yields G = 0.0364. With the above Tjj and G, from 
Eq.(0) we arrive at the recombination rate constant 

Orec = 1= a ~ 3.9 — a . (16) 

a/3 m m 

The dependence ^, instead of Orec oc a^, is a consequence of the recombination 

to a weakly bound s level and can be also obtained within the Jastrow approximation for 
the initial wavefunction: ipij = A/6xo(ri — r2)Xo(r2 — r3)Xo(i"3 "~ ri). This approximation 
was proved to be a good approach for atomic hydrogen and was later used for alkali 
atoms In our case, instead of Eq . (p!4D , we obtain tpij ^ V^Xoi^lxli^) and arrive 
at Eq.(|l6D, with 4 orders of magnitude larger numerical coefficient. Such a very large 
discrepancy occurs because both results are determined by distances x ~ a, where in our 
(rigorous) theory ipi is strongly reduced by the repulsive effective potential (see above). 
In the Jastrow approximation this reduction is not present. Moreover, ipij is resonantly 
enhanced at distances x < a. Thus, for large a the Jastrow approximation gives a wrong 
picture of three-body collisions and is absolutely inadequate to describe recombination to a 
weakly bound s level. 

The strong reduction of Orec due to the presence of a repulsive effective potential for 
three-body collisions can be treated as "quantum suppression" of three-body recombination 



(see related discussions in HT^J^l)- Nevertheless, a-^ec remains finite in the zero temperature 
limit. In fact, due to large values of a, it is rather large. It is also worth noting that for 
large and negative scattering length the quantity A(p) should have the form of a potential 
well, with a repulsive core at small p, and the picture of recombination collisions can be 
completely different. 

In trapped gases the kinetic energy of the third atom acquired in the recombination 
process usually exceeds the trap barrier, and such atoms escape from the trap. Thus, 
the loss rate for atoms is n = —Ln^, with L = Sctrec- For three-body recombination of 
ground-state He atoms Eq.(0) gives L ^ 2 ■ 10^^'' cm^/s. As the He-He interaction has 
Re ~ ISA^ a, this value of L is a very good approximation. More accurate calculation, 
using A(p) and Fx{p) determined for the He-He interaction (solid curves in Fig.2 and Fig.3), 
gives a correction of 10%. The same L is obtained for three-body recombination of spin- 
polarized He(2^S') atoms. In this case the result is less accurate, since the characteristic 
radius of interaction is somewhat larger (i?e ~ 35A). 

Qualitatively, the picture of an effective repulsion in three-body collisions, implying a 
strong reduction in the recombination rate constant, can be valid for systems with positive 
scattering length a ~ -Re- One can find such systems among the ultra-cold alkali atom gases. 
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is supported by the Royal Netherlands Academy of Arts and Sciences. This work was 
supported by the Dutch Foundation FOM, by NWO through project NWO-047-003.036, by 
INTAS and by the Russian Foundation for Basic Studies. 
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FIGURES 

FIG. 1. Three possible sets of coordinates for a three-body system. The relative coordinates 

are x, between two particles, and y, between their center of mass and the third particle. 

FIG. 2. The "effective potential" A as a function of pja. The solid curve is obtained from 
Eq.(8) using the ground state He-He potential (a = lOOA), and the dashed from Eq.(12). 

FIG. 3. The wavefunction F\{p/a) obtained from Eq.(9). The solid curve corresponds to A(p) 
for the ground state He-He potential, and the dashed curve to \{p) from Eq.(12). 
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